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The reaction pd — > 3 He 77 at threshold was used to provide a clean source of r\ mesons 
for decay studies with the WAS A detector at CELSIUS. The branching ratio of the 
decay 77 — > 7r + 7r~e + e~ is measured to be (4.3±1.3±0.4)xlCP 4 . 
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1 Introduction 

Radiative processes involving one or two photons are responsible for the most 
common decays of the lightest pseudoscalar mesons tt°, r\ and rj '. They are 
accompanied, as understood from Quantum Electrodynamics (QED), by a 
processes where a virtual photon converts into an electron-positron pair. The 
conversion decays are suppressed by a factor of the order of the fine structure 
constant a. The virtual photon probes the structure of the decaying meson 
and the interaction region in the time-like region of four-momentum transfer 
squared, q 2 , which is equal to the invariant mass squared of the lepton pair. 
An extensive review of the conversion decays of the light mesons is given by 
Landsberg [1]. 

Experimental information on the conversion decays of rj meson decays is scarce. 
The most frequent 77 decay into an e + e~ pair is 77 — ^ e + e _ 7 with a branching 
ratio (BR) of (6.0 ± 0.8) xlCT 3 [2]. Only a few hundred events were collected 
in recent experiments [3,4]. Another decay is ?7 — > 7r + 7r~e + e~ which is related 
to the radiative decay 77 -> tt+tt^ (BR = (4.69 ± 0.11) xl0~ 2 [2]). The first 
observation of one candidate event for this decay was reported from a hydrogen 
bubble chamber experiment by Grossman, Price and Crawford [5] in 1966. 
Recently semileptonic decays of the 77 mesons were studied by the CMD-2 
Collaboration [3] using the radiative decay <fi — > 777 as the source of rj mesons 
and 4 event candidates of 77 — > 7r + 7r~e + e~ were identified with an estimated 
background of 0.5 event. The average value of the BR extracted from the two 
experiments is (4.01! J) x !0~ 4 [2]. 

Predicted values for BR(rj — > 7r + 7r~e + e~) are given in at least three papers. 
Jarlskog and Pilkuhn [6] got an upper value of 3.1 x 10~ 4 . Faessler et al.[7] 
have calculated the decay rate within the Vector Meson Dominance (VMD) 
model and obtained a value (3.6 ± 0.6) xl0~ 4 . Picciotto and Richardson [8] 
got a value (3.2 ± 0.3) x 10~ 4 using a model which incorporates vector mesons 
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in the chiral perturbation theory Lagrangian. In this paper the influence of an 
intermediate p meson was taken into account both for pion and lepton pairs. 
The calculations are similar to the ones for the K L — > 7r + 7r~e + e~ decay [9,10]. 

Further interest in the 77 — > 7r + 7r~e + e~ decay comes from tests of CP violation 
in 77 decays which are often motivated by corresponding tests in Kl decays. 
A recent prediction (from CP violation) and observations of an asymmetry in 
the distribution of angles between the 7r + 7r~ and the e + e~ production planes 
in K L — > 7r + 7r~e + e~ [10,11,12] have triggered theoretical speculations that a 
similar observation in 77 — > 7r + 7r~e + e~ decay might reveal unexpected mech- 
anisms of CP violation in flavor conserving processes. It appears that there 
could be (hypothetical) CP violating contributions which are not constrained 
by the limits on the 77 — > tttt BR or the measured limit of neutron electric 
dipole moment. A measurement done with a sensitivity better than 10~ 2 for 
the asymmetry will provide a stringent constraint for such mechanisms [13]. 



2 The experimental method 

In the experiment performed by the CELSIUS/WAS A collaboration, the 77 
mesons were produced in the pd ^ 3 He 77 reaction at 893 MeV incident proton 
kinetic energy, very close to the threshold. This reaction was first employed 
as an 77 source for decay experiments at the Saturne II synchrotron at Saclay 
[14]. The production of 77 mesons was tagged by measuring 3 He at 0°. The 77 
production cross section in the reaction increases to a plateau value of 0.4 pb at 
about 2 MeV excess energy, where the background from prompt pd ^ 3 He7r + 7r~ 
reaction is at the percent level [15]. This tagging method enables to collect 
simultaneously an unbiased data sample of all 77 meson decays. 




Fig. 1. Cross section of the WAS A apparatus with Central Detector (CD) built 
around the interaction region, Forward Detector (FD) and Zero Degree spectrometer 
(ZD). The individual components are described in the text. 



In the CELSIUS/WASA experiment an internal deuterium pellet target was 
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used. The 3 He nuclei recoiling at 0° were detected in a zero-degree spectrom- 
eter (ZD) (Fig. 1). It used the accelerator quadrupole and dipole magnets to 
focus and deflect the 3 He nuclei onto a tagging telescope, which was placed 
inside the beam pipe 6.6 m downstream from the target [16]. The telescope, 
comprising two silicon detectors (total thickness 1.3 mm) and two germa- 
nium detectors (total thickness 29.4 mm), was able to stop 3 He nuclei up to 
400 MeV. Operated at liquid nitrogen temperature, the telescope provided an 
energy resolution of 1.5 MeV (FWHM) at 300 MeV [17]. The AE-E spectrum 
in Fig. 2 (left), where AE(E) is the total energy deposition in the two silicon 
(germanium) detectors, shows that the 3 He nuclei are easily separated already 
in the raw data. Fig. 2 (right) shows the energy spectrum of 3 He nuclei. The 
pronounced distribution between 280 and 305 MeV is due to r\ production. The 
two peaks, indicated by the structure at the top of the distribution, correspond 
to backward and forward emission of 3 He, respectively, in the center-of-mass 
frame. The smoother distribution is mainly due to direct 2-k production, yield- 
ing a background of less than 2% in the 7] meson tagging. The acceptance for 
3 He nuclei from the rj production is approximately 50% at 1 MeV above the 
threshold. 




E, MeV E, MeV 



Fig. 2. (Left) AE-E spectrum from the ZD, where AE(E) is the total energy depo- 
sition in the two silicon (germanium) detectors indicating very clean 3 He identifica- 
tion. (Right) Energy spectrum of 3 He nuclei. The pronounced distribution around 
300 MeV corresponds to the r? production. The cut-off at 200 and 380 MeV in the 
distribution associated with 2ir production is a result of the rigidity acceptance of 
the spectrometer. 

The rj decay products were detected in the central part of the WASA detector 
(Fig. 1) consisting of an electromagnetic calorimeter with over 1000 CsI(Na) 
crystals (SEC), a plastic scintillator barrel detector (PSB) and a drift chamber 
(MDC) made of thin-walled (25fim) aluminized mylar straw tubes placed in 
a 1 Tesla magnetic field produced by a very thin walled (0.18 X ) supercon- 
ducting solenoid [18] (SCS). The design of the WASA detector was optimized 
for measurements of decays with both electrons and photons. To minimize the 
probability for photon conversion it included a thin beryllium beam pipe with 
a wall thickness of 1.2 mm (0.0034 radiation lengths only). Fig. 3 shows an 
example of a candidate event for 77 — > 7r + 7r _ e + e _ decay registered in the CD. 
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Fig. 3. Example of a candidate event for rj — > 7r + 7r~e + e~ decay. Four charged tracks 
are measured in the MDC, all particles leave the volume of the chamber and deposit 
energy in the PSB and SEC (represented as the two outer rings). 

3 Analysis and results 

The candidate events for the 77 decay into 7r + 7r~e + e~ were selected by requiring 
four tracks caused by the charged particles in the central drift chamber with 
a possible energy deposit in the SEC and/or the PSB. Two particles should 
carry positive and two negative charge. A candidate event could in addition 
have a calorimeter hit-cluster not clearly associated with any of the tracks 
in the MDC (such a cluster will be referred to later in the paper as a "7"- 
cluster). These events could not be discarded directly without significant loss 
of acceptance (in the order of 35%). MC studies show that small but well 
separated additional clusters can arise from the 77 — > 7r + 7r _ e + e~ decay products 
due to secondary reactions in the calorimeter material resulting in emission of 
neutral particles. However, the energy deposit in such clusters is expected to 
be small. Only events with an energy deposit below 80 MeV in the "7" -clusters 
were accepted as candidates. 

Due to the clean tagging in pd — > 3 He 77, as discussed in the previous section, the 
only major background for the 77 — > 7r + 7r~e + e _ decay can be due to other 77 de- 
cays into two charged pions. A source of background could be the 77 — > n + n^n° 
decay with 7T° — > e + e~7 and the 77 — > 7r + 7r~7 decay with the photon undergoing 
conversion in the detector material. In order to reduce the latter background 
a vertex position constraint (Ax,Ay =±5 mm and Az =±50 mm) was used to 
ensure that all particles originated from the interaction region. That constraint 
removes more than 94% of 77 — > 7r + 7r~7 events with photon conversion. 

The two-lepton invariant mass distribution dT/dq 2 (i] — > 7r + 7r~e + e~) can be 
expressed approximately as dT/dq 2 = F(r) — > 7r + 7r~7) [QED] |F(g 2 )| 2 . The 
QED term ([QED]) is given by [1]: 
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where m e , are the electron and r] masses and is the invariant mass of 
the two-pion system. The q should be in the range 2m e < q < — 2m 7r ±. 
The QED term leads to a strong enhancement in the dT/dq 2 distribution at 
the lowest possible q values. The form-factor F(q 2 ) modifies the distribution 
mainly at large q. The most common assumption for the form-factor for r\ and 
7r° conversion decays is a dipole parametrization justified by p° dominance. 

The invariant mass is closely correlated to the opening angle between the 
leptons leading to a sharp peak at small opening angles. This feature of the 
process is used for particle identification. In order to distinguish between pions 
and electrons, the opening angle between two particles with opposite charges 
is calculated. The pair with the smallest opening angle is most likely the 
positron-electron pair. The efficiency of such a separation is demonstrated in 
Fig. 4 (left) using Monte Carlo (MC) simulation with the i] — > 7r + 7r~e + e~ 
mechanism given by the above formula. In addition an energy-momentum (E- 
P) method of particle identification is applied for all the particles which did 
not stop in the MDC. The energy deposit (dE) in the SEC and/or the PSB (E 
= dEsEC + dEps#) is combined with the momentum (P) calculated from the 
track curvature (the sign gives the charge of the particle). Fig. 4 (right) shows 
the simulated E-P distributions for all reconstructed tracks with associated 
energy deposition in the SEC. In that case the pion and electron bands, which 
are well separated, can be used to check and improve the identification based 
on the opening angle. The overall identification is correct in 90% of the cases 
for all tracks in the MC event. 




Z(7t + ,7T) ,(deg) P ,GeV/c 



Fig. 4. Particle identification for the MC sample of r\ — > 7r + 7r~e + e~ decay. (Left) 
Angle between the two pions (X-axis) versus the angle between the two leptons 
(Y-axis). The latter angle is smaller in 90% of cases. (Right) E-P distributions for 
all tracks with energy released in the SEC versus reconstructed particle momentum. 
The bands due to pions and electrons are well separated. 
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In order to suppress the background from the 77 — > 7r + 7r~7r° decays two con- 
straints optimized by the MC studies were used. The missing mass of the 
pions and leptons was constrained around the 3 He mass (in the range from 
2.6 GeV/c 2 to 2.85 GeV/c 2 ). This reduces the number of 77 — > 7r + 7r e + e and 
77 — > it + it~it events by 10% and 47% repectively according to the MC stud- 
ies. Some of the remaining events had so-called "7" -clusters in the calorimeter 
present (hit-clusters unassociated with a track in the MDC). Fig. 5 presents 
the invariant mass of the e + e~"7" system for the experimental sample with 
four charged particles and a neutral cluster in the final state. We find a clear 
peak around 7r° mass in accordance with the simulation of the i] — > n + n^n° 
decay (with n° — > e + e~7). Thus, for the class of events with a "7" -cluster, 
the invariant mass of the e + e~"7" system was calculated and an event was 
identified as a background event, if the mass was greater than 120 MeV/c 2 . 
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Fig. 5. The experimental invariant mass of e + e~"7" for the events with the 
"7"-cluster. The pronounced tt° peak in the experimental mass spectrum (black 
points) indicates a significant admixture of the background events. Superimposed 
are: a sum of the appropriately weighted MC simulation of rj — > ir + ir~e + e~, 
rj — ► 7T + 7r~7r° and 77 — > 7r + 7r~7 decays (full line) and rj — > 7r + 7r _ e + e~ decay (dashed 
line) . 

The invariant mass distribution, M n + n - e + e - , is plotted in Fig. 6 for the final 
data sample, first only for events without the "7"-clusters (left) and later us- 
ing all events (right). The experimental distribution is compared in both cases 
to the MC prediction taking into account the dominant background channels 
(77 — > 7T + 7r~7r° and 77 — > 7r + 7r~7) and the shape of the signal distribution. The 
contribution of the background events to the peak is estimated to be 4.6 and 
7.7 events (in the left and right plot respectively) by taking into account the 
known BR and acceptance from the MC studies. (The background contribu- 
tions are represented in the figures by the broken lines). The remaining events 
in the plots are attributed to the 77 — > 7r + 7r~e + e~ decay. This gives a number 
of reconstructed rj — > 7r + 7r~e + e~ events for the right plot in Fig. 6 equal to 
16.3±5.0 s t a t±2.0 syst (the systematic error includes the ambiguity of the back- 
ground contribution). The overall detection efficiency for rj — > 7r + 7r~e + e~ was 
estimated to be 16.5±0.2% where the error is dominated by the uncertainty in 
acceptance calculation due to extrapolation using different reaction models. 
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Fig. 6. The invariant mass distribution, M^+^-g+g-, plotted for events without the 
"7"-clusters unassociated with the MDC tracks (left) and for all events (right) (black 
points in both figures). The thick line represents the MC distribution corresponding 
to the sum of the signal (77 — > ir + ir~e + e~) and the background decays. In addition 
the contributions of the background from rj — > 7r + 7r~7r° and 77 — > 7r + 7r~7 decays are 
given by the dashed and dotted line respectively. 



The BR of 77 — > 7r + 7r e + e was obtained by normalizing to a sample of 77 — > n + n 
decays (5i? = (2.27 ± 0.04) x [2]) and 77 -> e + e"7 decays (5i? = (6.0 ± 
0.4) x 10~ 3 [2]) collected simultaneously. The data were analyzed using the same 
track finding processing and compatible cuts. This assures that many sources 
of errors, like uncertainties in particle identification and track finding effi- 
ciency, will be reduced or cancel. 



The selection criteria for the process 77 — > tt + 7t~7t° include requirement of two 
tracks with opposite charges identified as pions by the E-P method, two neu- 
tral clusters with energy E> 50 MeV each and an invariant mass of the two 
photons to be within the range 90 to 170 MeV/c 2 . Systematical uncertainties 
of the normalization were estimated by varying the constraints during sample 
selection. This lead to a change of the normalization factor within a 5% inter- 
val, what was included in the systematic error. The number of the identified 
events for this process was 17100±130 and the acceptance was estimated to 
be (33.4±1.9)% (the error is a cumulative systematic error of the normaliza- 
tion and includes the MC model dependence, the background and the sample 
selection) . 

A cross-check of the normalization with the 77 — > e + e~7 decay was done to test 
the MC predictions for the electron reconstruction efficiency. In the collected 
sample a clear peak of 270 events of the 77 — > e + e~7 decay has been observed 
after application of cuts on the charge balance, the overall missing mass and 
the invariant mass of the two leptons. Taking into account the calculated 
acceptance (20%) and assuming BR = (6.0 ± 0.8) x 10~ 3 [2] for the decay 
we got a normalization factor smaller by 5% from the one obtained with 
77 — > 7T + 7r~7r° decay, which is within the error of the given BR. The variation 
of cuts on the e + e~7 missing mass and/or replacing the cut on the lepton pair 



invariant mass by the cut on their relative angle lead to a variation of the 
normalization factor within 8% which is within one standard deviation of the 
number of i] — > e + e~7 events. 

Taking the obtained normalization factors and their statistical and system- 
atical uncertainties the result for the BR of the rj — > 7r + 7r~e + e~ decay is 
(4.3±1.3±0.4)xl(T 4 based on 16.3±4.9 stat ±2.0 S2/st identified events. 

It has been demonstrated that with the applied data selection the background 
channels are well understood. Our result increases notably the experimental 
information available for this r\ decay channel. The search for CP violation 
mechanism of [13], however, clearly requires much more data than what has 
been collected so far. 
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